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A sensitive search for CO J=1-0 emission in 4C41.17: 
high-excitation molecular gas at z=3.8 
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Abstract. We report sensitive imaging observations of the CO J= 1—0 line emission in the powerful high-redshift radio galaxy 
4C 41.17 (z = 3.8) with the NRAO Very Large Array (VLA ), conducted in order to detect the large concomitant H2 gas reservoir 
recently unveiled in this system bv lDe Breuck et ail 120051) via the emission of the high excitation J=4-3 line. Our observations 
fail to detect the J= 1—0 line but yield sensitive lower limits on the R43 = (4 - 3)/(l - 0) brightness temperature ratio of 
R43 ~ 0.55- S 1.0 for the bulk of the Ho gas mass. Such high ratios are typical of the high-excitation molecular gas phase 
"fueling" the star formation in local starbursts, but quite unlike these objects, much of the molecular gas in 4C 41.17 seems to be 
in such a state, and thus participating in the observed starburst episode. The widely observed and unique association of highly 
excited molecular gas with star forming sites allows CO line emission with large (high-J)/(low-J) intensity ratios to serve as an 
excellent "marker" of the spatial distribution of star formation in distant dust-obscured starbursts, unaffected by extinction. 
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1. Introduction 

The powerful radio galaxy 4C41.17 at z = 3.8 with a 
far-infrared (FIR) luminosity of Lfir ~ 10 13 L is the 
site o f an enormous starburst event in the di stant Universe 
(e.g. lEales et all 1 19931 iGrahametafl 1 1994. an interpre - 
tation put forth since its discovery JChambers etail [l990). 
This view received additi onal support from the discovery 
of a large dust mass JChini & Kriigell 1 19941 iDunlop et all 
11994 . and its concomitant molec ular gas reservoir with 
M(H 2 ) ~ 5 x 10 10 M dPe Breuck et all I2005I) . which 
presumably fuels the star formation. This follows the path 
of similar discoveries made in the last decade, uncovering 
vast molecular gas masses residing in a variety of high- 
redshift objects. These include QSO s dBarvainis et all 1200' 



Omontetal., 1996; Ohtaetal., 1996), submm-selected galax- 
ies IFraver et all 11 998: Ner i et all 120031) . Ly-b reak galaxies 
( [Baker et all l2004 and powerful radio galaxies IKlamer et all 
(2005). The highest redshift at which CO has been detected 
stan ds now at z = 6.42, when the Universe w as less than a Gyr 
old dWalter et all 120031 iBertoldi et alll2003l) . 

Aside from the intrinsic interest of studying the molecu- 
lar gas in distant powerful radio galaxies, the likely formation 
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sites of pr esent-day giant elliptica l s and their surrounding clus- 
ters (e.g. IPentericci etail 1 19971 lArchibald etail 120011) . the 
case of 4C41.17 is special since it was the first non-lensed 
high-redshift object for which dedicated, sensitive observations 
were c onducted to detect its m o lecula r gas mass Jlvison et all 
119961 E arvainis & Antonucci, 1996). Observations of CO 
J=l-0 are of particular importance since they provide the 
least excitation-biased method of tracing metal-rich H2 gas 
(Ei/kfi ~ 5.5 K, n cl ~ 400 cm -3 ), while most of the reported 
searches/detections of molecular gas in high-redshift objects 
use C O J+l— >J, J+ 1 >3 and usually J _ 4-3 or higher transitions 
(e.g. lEvans et ail Il996l Ivan Oiik et all ll997tlDownes et ail 
1999). These transitions, with their much higher critical den- 
sities and energy levels (e.g. for CO 4-3: n cl ~ 2 x 10 4 cm -3 , 
E4/kfi ~ 55 K) trace much warmer and denser gas, and thus 
may not reveal the distribution of the bulk of the molecular gas. 



Evidence for molecular gas, sub-thermally excited over 
large scales, has emerged for a few objects at high r edshifts 
(e.g. IPaoadopoulos & Ivisoiil2002llGreve et alll2003l) and for 
this reason we followed up the recent detection of CO J = 4-3 
in 4C41.17 with sensitive observations of the J = 1-0 transi- 
tion. At z = 3.8 the latter is redshifted to 24 GHz and becomes 
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accessible to the K-band receivers of the VLA. 1 Our sensitive 
observations, which unlike past efforts had the advantage of 
knowing the line center from the detection of CO J = 4-3, do 
not detect any J = 1-0 emission but the lower limits they set on 
the R43 = (4 — 3)/(l - 0) line ratio imply a high excitation for 
the bulk of the molecular gas. This suggests that, unlike local 
counterparts, most of the molecular gas reservoir in this distant 
galaxy participates in the observed starburst event. Throughout 
this work we assume H = 75 km s _1 Mpc~' and q D = 1/2. 



2. Observations 

The observations of 4C41.17 were conducted with the VLA 
in its D configuration with four ~8-hr tracks obtained during 
2004 August 12-22 UT and a phase-tracking center at: RA 
(J2000) = 06 h 50 m 52U50 and Dec (J2000) = 41° 30' 30".801. 
The CO J = 1-0 line (v rest = 115.2712GHz) at the redshift 
Z = 3.7978 of the radio gala xy (centered on its CO 4-3 emis- 
sion; |B^Briu£kii^il H 005 ) shifts to v obs = 24.0258GHz, ac- 
cessible to the VLA's K-band receivers. Aiming solely at the 
detection of line emission rather than obtaining information 
about the velocity distribution of the molecular gas, we used the 
continuum correlator mode with one IF pair (left and right cir- 
cular polarization) tuned to observe the line+continuum while 
the other one simultaneously observed the adjacent continuum. 
At 24 GHz the effective IF bandwidth, Av IF = 45 MHz, cor- 
responds to AVif ~ 560 km s . The wide velocity range 
(~ 1000 km s" 1 ) revealed by the CO J = 4-3 observations en- 
sures no sensitivity loss due to velocity dilution of the J = 1- 
line within the IF bands. Two separate tunings were neces- 
sary (and were equally distributed per track) to observe the 
entire expected line width. These [on-line, continuum] tun- 
ings, optimized to observe the 'blue' and the 'red' veloc- 
ity components of the CO J = 4-3 emission, were: [24.0649, 
24. 1649] GHz and [24.0351, 23 .9351] GHz, corresponding to 
[-490, -1735] km s -1 and [-115, +1130]kms- 1 velocity off- 
sets relative to v Q b s (z). 

During our observing sessions the system temperatures 
were 7\,v S ~ (70 - 90) K, and a fast-switching technique was 
used JCarilli & Holdawavt [l999), recording data every 5 sec 
(integration time per visibility) with 20 sec spent on the com- 
pact calibrator 06465 +445 1 3 and 150 sec on source; the typical 
slewing time to the calibrator, 3.5° away, was 20 sec. This en- 
sured very good amplitude and phase calibration with a result- 
ing phase |0 rms | 55 25° for most of the time. The latter was esti- 
mated from the residual phases of our amplitude/phase calibra- 
tor after applying antenna gain solutions smoothed to ~ 3xAt ca i 
(At ca i ~ 3.2 min: interval between successive observations of 
the calibrator). The flux density scale was set using 3C286 
(2.59 Jy) and 0713+438 (0.55 Jy). The resulting bootstrapped 
flux of 06465+44513, used for the final amplitude/phase cali- 
bration of the data, was S (24 GHz) = 3.15 Jy (averaged over 
all four days), and carries an uncertainty of ~ 15%. 
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Fig. 1. Radio continuum map of 4C41.17 at 24 GHz re- 
sulting from imaging all the uv datasets (see text). Robust 
weighting was applied (with robust = 0) and clean was 
employed, converging to a total flux Scl = 4.88 mJy). 
The map was then convolved to circular synthesized beam 
of FWHM=2.5 shown at the bottom right. The contours 
are at (-3,3,6,9,12,15,18,30,50,70,80) x cr, with cr = 
28 yuJy beam 1 . The cross marks the position of the radio 
core, RA (J2000): 06 50 52.15, Dec (J2000): +41 30 30.80 
llCarillietalJ. fi 994). which is also the phase tracking center of 
the observations presented. 



2.1. Data reduction: the search for CO J = 1-0 line 
emission 

The measured visibilities were carefully edited and calibrated 
using standard techniques described in the NRAO ATPS 
Cookbook and then Fourier transformed (using the task imagr) 
to produce images. CLEANing was applied to produce the final 
(deconvolved) images of the emission of 4C41.17 at 24 GHz. 
In all our naturally- weighted maps the measured r.m.s. noise 
was close to the theoretical value expected from the number 
of the visibilities imaged, the system temperatures during our 
observations, two IFs per frequency and an assumed antenna 
aperture efficiency of r\ a = 0.40 at 24 GHz 2 . A map of the 
entire field of view of the primary beam (~ 2 ), using all 
the data is shown in Fig. ^ The flux densities measured for 
the core and its SW feature are: S core = (2.97 ± 0.45) mJy 
and Ssw = (0.73 + 0.12) mJy, in good agreemen t with those 
reported for these components by llvison et"al] ill 99 6). The 
SW feature corresponds to feature 'A in lo w-frequency maps 
(Chambe rs et all Il990t ICarilli et all fl994) and is dominated 
by optically thin synchrotron emission from th e radio lobe, 
with a steep spectral index a = -1.7 (S v oc V) Jlvison et all 
1996). For the emission feature at AO 50 to the South (also 
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06 50 52.6 52.4 52.2 52.0 51.8 51.6 51.4 51.2 

RIGHT ASCENSION (J2000) 

Fig. 2. Radio continuum map of 4C41.17 at 24 GHz (con- 
tours) overlayed on the 2.15-/im image obtained at the Keck 
Observatory (Graham et al. 1994). The contour levels are (3, 
4.25, 6, 8.5, 12, 17, 24, 34, 48, 68, 96) xcr, with cr=26 fiiy 
beam -1 . 



seen in the maps of Carilli et al. 1994) we find a probability of 
it due to chance of P(A0) ~ 0.25, thus the association is not 
significant. 3 

An overlay of our new 2 4-GHz map with those obtained at 
1 .4 and 5 GHz JCarilli et allll994l) reveals excellent correspon- 
dence of the emission features across all three wavelengths, a 
testimony to the good phase calibration at 24 GHz resulting 
from our fast-switching scheme. A deep K$ -band (2.15//m) 
image shows a faint feature just upstream from the SW radio 
lobe A (see Fig. |2J. At a redshift of z = 3.8, 2.15^/m cor- 
responds to rest-frame Z?-band, expected to be dominated by 
young stars. This is especially true since the K$ -band is free 
of emission-line contamination, and the scattered AGN light is 
expected to be negligible ~35 kpc from the radio core position. 
If this component is indeed at the redshift of 4C41.17, it would 
indicate a second locus of jet-induced star formation in th is dis- 
tant radio galaxy JPev et allll997tlBicknell et al.Ll2000h . 

We searched for CO J = 1-0 emission by comparing the im- 
ages "on" and "off" the expected line for the two different fre- 
quency settings covering its velocity range (as revealed by the 
CO J = 4-3 emission). We did not find any significant excess 
flux density in the maps "on" with respect to those "off" the 
line, within the noise levels achieved, while the total CLEANed 
flux density of the radio galaxy is similar for these maps and 
for both frequency settings. We produced line-only maps us- 
ing the AITS task uvsub: the clean components of the con- 
tinuum maps were Fourier transformed back to the visibility 
plane, re-sampled by the uv sampling function (identical for 
the simultaneously obtained on-line and off-line datasets), then 



3 We estimate P (A9) following the method outlined by 
iDownes et al] jl9Bfit) and the 1.4-GHz number counts reported 
bv lHopkinsI <2003l) after we extrapolate a flux of S l 4 ghz - 5 mjy for 
our source (assuming a spectral index of a = -0.7 between 1.4 and 
24 GHz). 



subtracted from the (line + continuum) visibilities. The result- 
ing continuum-free visibilities were then Fourier transformed 
(without CLEANing) to produce naturally weighted (for max- 
imum sensitivity) images (Figs[3]and0J. Significantly tapered 
images from these visibilities were also made to investigate the 
presence of any extended, low-brightness, line emission but no 
evidence of this was found. 

Tantalizingly, in the frequency setting covering the "blue" 
profile of the CO J = 4-3 line some faint emission can be seen 
in the line map, located between the two radio continuum com- 
ponents (Fig. However, this is not statistically significant, 
especially since the noise measured within the area containing 
the radio galaxy (and where the clean components of its con- 
tinuum are derived) is somewhat higher than in the rest of the 
image, namely cr eff = 65 /iJy beam -1 . This is expected since 
subtracting the noisy clean components, Ncl, defining the ra- 
dio continuum emission will add noise to that part of the re- 
sulting continuum-free image (<x eff < V2cr ~ 80 ji]y beam -1 , 
simply because Ncl < iVpbteis). We use cr e ff to express the 1-cr 
upper limit of the velocity-integrated CO J = 1-0 line from 



AS, _o = f (5S v dV = c f^E) VM>e ff , 

JAv \ Vobs I 



(1) 



where Nb is the number of beams associated with the puta- 
tive CO-emitting area. A single IF covers most of each of the 
detected J = 4-3 "blue" and "red" velocity components, with 
AS i_o = 0.035 Jy km s -1 (point source sensitivity, Nb = 1), and 
a corresponding limit on the J=l-0 luminosity per component 
of L co ~ 10 1() Kkms -1 pc 2 . The corresponding lower limit for 
the velocity/area-averaged R43 line ratios can be found from 
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The CO J = 4-3 emitting regions are unresolved by the 
#hpbw ~ 6" of the Plateau de Bure Interferometer, thus 
Nb(max) ~ 3.5. Using the numbers reported by De Breuck et 
al. for $4_3, and considering a 2-cr upper limit for the J = 1-0 
emission, with N b = 3.5, yields > 0.30 and R® > 0.55 
for the "blue" and "red" component, the latter containing most 
of the H2 gas mass in this system. For a CO emitting region 
smaller than our resolution limit of ~ 3", Nb = 1 and these lim- 
its become R^ 0.55 and ^ 1 ■ These are the highest min- 
imum values compatible with the sensitivity of our measure- 
ments, and the most like ly ones since in this object f? co ^ 0.3" 
(De Breu ck et all 12005). and are adopted in all our subsequent 
analysis. 



3. Discussion 

The models of the molecular gas physical condit ions con- 
strain ed by single or even multiple line ratios (e.g. lAalto et all 
1 1995b contain significant degeneracies, but for modest-to-high 
values of CO line ratios with large J-level differences (and thus 
difference between n cl and E uppel /ke) these are significantly re- 
duced. The CO J = 4-3 transition in particular marks a broad 
excitation turnover beyond which line flux densities and thus 
the line ratios Rj+u (J + 1 > 4) become highly dependent on 



4 



Papadopoulos et al.: CO J=l-0 in 4C41.17 



the ambient gas conditions and can vary widely (Nieten et al., 
1999; Carilli et al., 2002). Moreover, because of the difficulty 
of routine observations at v S V43 ~ 461 GHz through the 
Earth's atmosphere, ratios like R43 are sensitive probes of gas 
excitation that are rarely available in the local Universe. 

The lower limit of R43 > 0.55 of the "blue" emission 
component in 4C 41.17 is within the range of the R43 val- 
ues measured in another powerful radio galaxy, 4C 60.07 at 
z~3.8 JGreve et all l2004h . and suggests highly excited gas. 
By comparison, the globaly averaged CO emission for the 
Galaxy measured by COBE, yields R 43 ~ 0.10 - .17, indica- 
tive of a quiescent gas phase dFixsen et all Il999l) . A single- 
phase Large Velocity Gradient (LVG) model, restricted by set- 
ting T k = T dust ~ (50 - 60) K (De Breuck et al., 2005), yields 
n(H2) <; 10 3 cm -3 for most solutions, which is on the upper end 
of the mean densities of Giant Molecular Clouds in the Galactic 
disk (~ (10 2 — 10 3 ) cm -3 ). If an ensemble of virialized clouds is 
assumed, the average cloud density and its associated velocity 
gradient are no longer independent but related as 

where a = 0.5 - 2.5 depends primarily on the cloud den- 
sity profile (Papadop oulos & Seaauisl Il999l) . For dV/dR = 
(dV/dR) VIR (and an abundance [ 12 CO/H 2 ] ~ 10~ 4 ), the LVG 
solutions compatible with R43 ^ 0.55 have ~ (40 - 50) K 
for densities of ~ (1 - 3) X 10 3 cm" 3 . However, in starburst 
environments highly non-y irial motions have been deduced 
llDowne S &Solomonlll998h and thus LVG solution constraints 
based on known dust temperatures are better than assumptions 
of virialized I2 CO-emitting clouds. For the "red" kinematic 
component, with R43 ^ 1, our models yield n(H 2 ) ^ 10 4 cm -3 
(for Tk ~ (50 - 60) K), typical for H 2 gas in star-forming 
regions. The solutions compatible with R43 > 1 and T^ ~ 
(50 - 60) K typicaly have moderate CO J=l-0 optical depths 
(tio ~ 1 - 2), a condition encountered often for t he warm and 
turbu lent molecular gas in starburst environments jAalto et all 
119951) . 

3. 1 . The dominance of high-excitation gas in 4C 41. 17 

One-phase LVG modeling of CO line ratios encompassing en- 
tire galaxies (the case for most high redshift CO line detec- 
tions), can yield only indicative results regarding which con- 
ditions are prevelant in their large molecular gas reservoirs. A 
different way to determine the relative prevalence of a highly- 
excited, starburst-fueling, gas phase is to make use of its well- 
established characteristics (and typical CO line ratios) along 
with those of a quiescent, star-formation-idle one. These two 
phases have been extensively studied in the local Universe 
where observations of several molecular and atomic lines, as 
well as a high spatial resolution, allow their detailed and realis- 
tic modeling (see Papadopoulos, Thi, & Viti 2004 for a review). 

In the local Universe, the unique association of Rj + i j ~ 1, 
J + 1 < 4 ratios (typical of high-excitation molecular gas) 
with sites of vigorous star formation has become apparent 
since the first extragalactic detection of CO J=4-3 emission 



by iGiisten et all l[l996) towards the star-forming regions of 
the nearby archetypal starburst galaxies M82, NGC253 and 
IC342 where typical (4-3)/(2-l) CO line ratios of R 4 3 /2 i £ 
0.5 - 0.8 were measured. For M82 this wa s confirmed by 
much more detailed studies llMao etall EoOOi. while a recent 
study demonstrated the existence of an additional massive low- 
excitation gas phase extending much further (~ 3 k pc) than the 
starbu rst-related gas residing in its inner 500 pc JWeifi et all 
1200.4 . 

For starbursts at high redshifts, the distribution of star form- 
ing regions is one of the crucial aspects that can distinguish 
between scaled-up (in terms of star formation rate) versions 
of local starburst galaxies (usually with compact star-forming 
regions) and very different, galaxy-wide events. Extended CO 
line emission with large (high-J)/(low-J) intensity ratios can be 
an excellent "marker" of the spatial distribution of star forma- 
tion in distant dust-obscured starbursts, unaffected by extinc- 
tion. In the case of 4C41.17 the CO J=4-3 emission is unre- 
solved, but the distribution of a high-excitation, star-forming, 
gas phase and its concomitant starburst relative to the total CO- 
emitting area can be constrained from our limit on the global 
R43 ratio. Indeed, for the large velocity gradients characteriz- 
ing molecular clouds (thus no significant radiative "shadow- 
ing" between different gas phases) it is 



_ R^+PhAiR^ 



where R£ J and are the (4-3)/(l-0) ratios of the 

velocity/area-averaged, CO emission of the high-excitation and 
low-excitation molecular gas phase, fhi = i^h/^i the ratio of 
their emitting areas, and p hl = [J(T[ h) )-J(T bg )]/[J(Tf)-J(T bg )] 
(for optically thick and thermalized J=l-0 transition for both 
phases), where J(T) = hv/k B (e hv/kBTk - l)" 1 (v = vio ~ 
115 GHz) and T bg = (1 + z)T cmb ~ 13K is the CMB tempera- 
ture at z ~ 3.8. Setting R^ = 1 (typical for the high-excitation, 
star-forming gas phase e.g. WeiB et al., 2005 and references 
therein), and R^ = 0.15 fo r the quiescent and star-forming- 
idle one dFixsen et all Il999l) . the lower limit R43 ~ 0.55 cor- 
responds to phifhi ~ 0.9. For typical temperatures usually asso- 
ciated with these gas phases: T® ~ 20 K and T^. h) ~ 60 K, it is 
f hi ~ 0.13. The latter amounts to f h = f h i/(l + fhi) ~ 0.10 of the 
total CO-emitting region occupied by the high-excitation gas 
phase and its concomitant starburst (by comparison in M 82: 
fh ~ 0.03). A maximum value for fhi and fh can be found by set- 
ting T^ h) ~ t£* (thus phi ~ 1), i.e. the difference in the R43 ratio 
(and the excitation state of J=4-3) between the two gas phases 
is solely due to a difference in densities, yielding fhi ~ 0.9 and 
fh ~ 0.47. In practice both temperature and density effects play 
a role and the true fh will lie between the aforementioned es- 
timated values. These are applicable to the "blue" component; 
in the "red" component, which hosts the radio-loud AGN, the 
R43 ^ 1 limit implies a high-excitation gas phase that fully 
dominates the CO emission. 
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4. Conclusions 

We report sensitive new observations of the CO J = 1-0 tran- 
sition in the distant, powerful radio galaxy 4C41.17 at z = 3.8 
using the NRAO Very Large Array. Our results can be summa- 
rized as follows: 

1. The CO J=l-0 line emission remains undetected but the 
new sensitive upper limit and the detected J =4-3 line yields 
a brightness temperature ratio of R43 ~ 0.55- ^ 1 for the 
bulk of its H2 gas reservoir. This ratio is ~ 3-6 times higher 
than the one in quiescent environments of the Galaxy, and 
is typical of starburst environments. 

2. Single-phase Large Velocity Gradient models, under the re- 
striction of T k = T dust , yield n(H 2 ) ~ (10 3 - > 10 4 ) cnT 3 , 
more typical of the high-excitation gas "fueling" local star- 
bursts rather than a more diffuse, low-excitation, phase. 

3. Assuming the molecular gas reservoir to be segregated be- 
tween a high-excitation star-forming, and a quiescent star- 
forming-idle phase, with R43 ratios and gas temperatures 
typical for the local Universe, we find the former to oc- 
cupy at least ~ 10% and up to ~ 100% of the total CO 
emitting area. This is much higher than what is found in 
local starbursts like M82, and equivalently demonstrates 
the predominance of a high-excitation gas phase fueling an 
extreme star formation event in this distant powerful radio 
galaxy. 
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Fig. 3. Naturally-weighted maps of ("blue" line)+continuum (top left), continuum (top right), and line only (bottom) emission 
(see text). Contours: (-3, -2, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 25, 35, 45, 55, 65, 75, 85) x a (cr = 55/Jy beam -1 ). The restoring beam is 
3.27" x 3.16" and is shown at the bottom right corner (clean was not applied on the continuum-free map). 
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Fig. 4. Naturally-weighted maps of ("red" line)+continuum (top left), continuum (top right), and line only (bottom) emission (see 
text). Contours: (-3, -2, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 15, 25, 35, 45, 55, 65, 75, 85) x cr (a = 55yuJy beam -1 ). The restoring beam is 
3.3 1 x 3.13 and is shown at the bottom right corner (clean was not applied on the continuum-free map). 



